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The stationary isotropic Poisson line n etwork was used to derive upper bounds on 
mean excess network- geodesic length in lAldous and Kendalll ( 20081 ). This new paper 



presents a study of the geometry and fluctuations of near-geodesics in such a network. 
The notion of a "Poissonian city" is introduced, in which connections between pairs of 
nodes are made using simple "no-overshoot" paths based on the Poisson line process. 
Asymptotics for geometric features and random variation in length are computed for 
such near-geodesic paths; it is shown that they traverse the network with an order 
of efficiency comparable to that of true network geodesies. Mean characteristics and 
limiting behaviour at the centre are computed for a natural network flow. Comparisons 
are drawn with similar network flows in a city based on a comparable rectilinear grid. 
A concluding section discusses several open problems. 
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1 Introduction 



The " Poissonian city" is a network of connections based on a Poisson line process. lAldous and Kendal] 



( 120081 ) used such a network to address a problem in frustrated optimization: construct planar 



networks connecting a large number of nodes such that 

1. the total connection length is not much larger than the minimum possible connection 
length, but also such that 

2. the average connection distance between two randomly chosen nodes is not greatly in 
excess of the Euclidean distance. 

It transpires that networks satisfying criterion [1] may be augmented by sparse Poisson line 
processes so as to satisfy criterion [2j More precisely, suppose that n nodes are distributed in 
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an arbitrary fashion (deterministically or randomly) over a square of total area n. Recall that 



(Promel and Steeer 


1997, 


Yukich 


. 1998) 



Aldous and Kendall! (120081 . Theorem 1 (b)) show that augmentation by 



a sparse Poisson line process can convert a Steiner minimum tree into a network whose total 
connection length is only slightly increased, but which now delivers a mean connection dis- 
tance which is no more than O(logn) in excess of Euclidean distance. U nder a suitable weak 
unifor mity condition on the empirical spatial distribution of the nodes, lAldous and Kendall 
(120081 . Theorem 2) also establish a lower bound on the mean excess; it must be of order at 
least fl ( \/Tog n). 

The primary motivation of this previous work was to gain better understanding of the 
behaviour of network statistics (such as the mean excess network length) for entirely general 
networks. However the appearance of Poisson line processes in the upper bound result moti- 
vates a more detailed study of the "Poissonian city" generated by a unit intensity stationary 
isotropic Poisson line process. What can be said about the "near-geodesics" used to establish 
the upper bound? how close are they to true geodesies? how does random fluctuation affect 
excess length? and what about traffic flow on such a network? These questions are addressed 
below; their answers require the use of Levy subordinators and self-similar Markov processes 
(something of a novelty in stochastic geometry), and a curious improper anisotropic Poisson 
line process. 

Previous relevant work includes: the note by lDavidsonl (119741 ). who gi yes a qualitative ar- 
gumen t showing that the Poisson line process provides good connections; iRenyi and Sulanke 
( 119681 . Satz 5), who derive a result similar to the mean-excess result but concerning numbers 
of edges rather than length, and based on a fixed number of random lines; and recent h igher- 
dimensional generalizations of the Renyi -Sulanke work by iBoroczkv and Schneider! (12008 . 
Theorem 1.3). We also mention work by IVoss. Gloaguen. and Schmidtl (120091 ) on limit dis- 
tributions of shortest paths from subsidiary to major nodes in hierarchical networks based on 

rando m tessellations. Finally we note the interesting work of lBaccelli. Tchoumatchenko. and Zuyev 
( 120001 ). related to the concept of spanners from graph theory (a geometric spanner is a 
planar graph connecting a set of nodes for which the graph distance between any two 
poin ts is less than some fixed m ultiple of Euclidean distance; see for example the e xposi - 
tion 



Narasimhan and Smidl . 120071 ) . The networks constructed in lAldous and Kendalll (120081 ) 



are averaged rather than uniform versions of geometric spanner networks, for which the fixed 
multiple of Euclidean distance is replaced by a logarithmic additive excess and a specific 
constraint is imposed on the total network length (rather than, say, small vertex degree or 
total number of edges). 

In the remainder of this introductory section we introduce basic notation and concepts, 
and enumerate the questions which will be addressed concerning the behaviour of near- 
geodesics and traffic flow in the Poissonian city. 



1.1 Notation and basic concepts 

We begin by presenting a brief summary of stationary isotropic Poisson line processes, so 
as to fix notation and to collect some facts about line processes which will be used below. 
Further information can be found for example in IStoyan. Kendall, and Meckel (119951 ). The 
ensemble of undirected lines i in the plane may be viewed as a Mobius strip of infinite width, 
or as a once-punctured projective plane (since such lines can be produced as intersections of 
planes through the origin in 3-space with the plane z — 1, in which case the plane through 
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the origin and parallel to z = 1 does not produce an intersection). It is often convenient to 
parametrize this ensemble of lines i by representing lines using points (r, 9) where r is the 
perpendicular signed distance from the line I to a reference point, and 9 6 [0, ir) is the angle 
that i makes with a reference line running through the reference point. A unit intensity 
stationary isotropic Poisson line process ( "Poisson line process" for short) is determined as 
a Poisson point process on the representing projective plane using the intensity measure 
| dr d9. The factor | ensures that the mean number of Poisson lines hitting a line segment 
is equal to the length of the segment. 

Slivynak's theorem on the Palm distribution of a Poisson process applies here: if we 
condition on a specific line I belonging to the Poisson line process then the residual line 
process is still unit intensity stationary isotropic Poisson. 

An alternative parametrization (x, 9), sometimes of use, employs the line angle 9 as above 
and x the signed distance along the reference line from the reference point to the intersection 
of i with the reference line. This representation breaks down when 9 = (not a major issue in 
the case of a Poisson line process, for which the set of lines at 9 = has zero probability). In 
these coordinates the unit intensity measure is ^sm9dx d9; the sine- weighting corresponds 
to a length-biasing phenomenon when sampling Poisson lines according to their intersections 
with a test line. In particular, if two lines are conditioned to pass through a given point 
then they form an exchangeable pair: one having uniform direction, and the angle a e [0, 7r) 
between them having density |sina independent of the direction of the first. 

Viewed as a random measure, the Poisson line process generates a measure via the mean 
total length of lines intersected with a given set. Testing against a unit disc, we can compute 
the resulting length intensity as |. It follows from the above (and in particular from Slivynak's 
theorem) th at the point process of intersections of lines from the Poisson line process has 
intensity | (IMilesl . Il964j . Theorem 2). 




source p 



destination p + 



Figure 1: A caricature of the procedure of finding a route using a Poisson line process; consider 
only those routes involving the two vertical lines I and one of the other lines. Here a possible 
route is indicated by arrows. 



The following caricature supplies a good intuition as to where the logarithmic excess might 
be located on a typical route on a network based on a Poisson line process. Consider a network 
formed between just 2 nodes p~ (the source) and p + (the destination), with lines provided 
by a unit rate stationary isotropic Poisson line process IT. Let the 2 nodes be separated by 
distance n. We condition the line process to contain two lines I running through source and 
destination nodes which are both perpendicular to the segment connecting p~ to p + (Figure 
[TJ. We consider only those routes which involve the conditioned lines £~ , respectively £ + , 
running through p~, respectively p + , and also just one other line of the Poisson line process. 
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Consider the set of lines £ which intersect both l~ at distance at most Cy/n from p~, and 
also £ + at distance at most Cy/n from p + . Classic stochastic geometry arguments (based on 
inclusion-exclus i on an d a special case of the "Buffo n-Sylvester problem" — see for example 
Ambartzumianl . Il990l ) then show that the invariant measure of this line set is given by half 
the difference between the summed length of the two green lines minus the summed length 
of the two red lines in Figure EJ 



2 c Vn 





n 













2c Vn 



Figure 2: Illustration of a classic stochastic geometry construction for calculating the invariant 
measure of the set of lines hitting both of the two vertical line segments £^ marked out by the 
two horizontal lines. 



Hence the probability of no unconditioned Poisson lines falling in this set is 
exp (-§ (2\/4c 2 n + n 2 - 2n\} > exp (-2c 2 ) , 
and as a consequence the resulting mean excess is bounded below by 



i~ 2c2 dc 



Tin 




attributable to the parts of the route which lie on the conditioned lines l~ , £ + . (Excess 
along the unconditioned line £ itself is bounded above by \/4c 2 n + n 2 — n < 2c 2 and hence is 
negligible in the case of large n.) 

This rather trivial example makes it clear that being permitted to use more than one line 
(in addition to the two conditioned lines) will reduce the excess principally by rounding off 
the corners at st art and finish of the j o urney . Thus it is clear (as is indeed apparent from the 



computations in lAldous and Kendall! . 120081 . Theorem 3) that the logarithmic excess in the 
full construction is a cost which arises entirely from the business of getting onto and off an 
efficient route between source and destination. 



1.2 Making connections 

Any two specified points p~ and p + in the plane will almost surely not be hit by any of 
the lines of a given isotropic stationary Poisson line process IT, and therefore will fail to be 
connected by II. Accordingly we establish the convention that movement from p~ to p + 
occurs as follows: first use the Poisson tessellation to construct the cell C(p~ ,p + ) containing 
p~ and p + which arises by deleting all Poisson lines which separate p~ from p + . Now proceed 
from the source p~ in exactly the opposite direction to that of p + , till one first encounters a 
Poisson line (which will be part of the cell boundary dC(p~ ,p + )). Then continue along the 
line in one or the other direction, clockwise or ant i- clockwise, proceeding along the boundary 
of the cell C(p~,p + ). Continue until one reaches the ray extending from p~ and through 
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p + . Then proceed down this ray to the destination p + (Figure [3]). Thus we consider near- 
geodesies; routes based on semi-perimeters of the cell C{p~ ,p + ). 

Evidently this is a conservative option for plumbing nodes in to the Poisson network 
produced by II, suitable if we wish to produce upper bounds on connection lengths and 
adopted without further comment in what follows. We suppose the choice of whether to travel 
clockwise or ant i- clockwise round the cell C(p~,p + ), (equivalently, which semi-perimeter to 
choose) is made at random and equiprobably, independently for each pair of nodes p~ , p + . 
(As mentioned in Section interesting and hard problems arise if the choice of route for a 
specific pair is influenced by the flow in the entire network.) 

In contrast to true geodesic connections, these routes can be viewed as outputs from 
an unsophisticated but direct semi-perimeter algorithm: if one is on a Poisson line and 
encounters another Poisson line then one chooses (from the three onward paths) that path 
which leads closest to the eventual destination without overshoot. This focuses attention on 
the Poissonian city, a region connected by routes based on a fixed stationary isotropic Poisson 
line process and following the above convention so as to ensure that the line process actually 
connects nodes. Que stions addresse d in Section [2] of this paper, filling in and extending the 



results announced in iKendalll (120081 ). include: 



1. What can one say about the basic geometry of these routes? Computations from 
Aldous and Kendall! (120081 . Theorem 3) yield | logdist(p~,p + ) as asymptotic mean ex- 



cess length as distQo ,p + ) tends to i nfinity. This can be viewed as a quantitative 



development of the announcement by iDavidsonl (11974 . Theorem 5(ii)); but by how 



much does the travelled path deviate laterally from the Euclidean connection, and at 
what point is that lateral deviation greatest? (See Theorem [TJ) 

2. What is the order of random variation of the route-lengths? (See Theorem [31) 

3. What might be said about how actual geodesies differ from these routes? (This is 
discussed in Section 1231 ) 

4. Finally, might actual geodesies produce substantially smaller mean excess lengths? 
(Theorem H] shows that the mean excess of true geodesies is comparable to the mean 
excess of semi-perimeter paths.) 

1.3 Traffic flow 

Given a Poissonian city, it is natural to consider traffic flow. Suppose for example that the 
city is represented by ball(o,n), a disk centred at o and of radius n, connected by roads 
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provided by a stationary isotropic Poisson line process II. Suppose that each pair of points 
p~ and p + in the disk generates a constant infinitesimal amount of traffic dp~ dp + divided 
equally between each of the two connecting routes generated according to the semi-perimeter 
algorithm described above. Suppose further that we condition on the event of a Poisson line 
passing through the centre o of the disk. Let 

V n = {(p~,p + ) Gball(o,n) 2 : p\ < p+ , oedC(p~,p + )} (1) 

denote the region in 4-space corresponding to p~, p + in ball(o, n) for which p~ lies to the left 
of p + (imposed by the inequality p~{ < pf, where p~ , p + are the ^-coordinates of source and 
destination nodes) and one of the two possible routes passes through o. 
Questions addressed in Section [3] include: 

1. What is the dependence on n of the mean E [T n ] of 

Tn=1 2 l l(p-,P + )tv n ] dp' dp + = | / / l [ Pl<P +, oedC( P -,p+)} dp~dp + , 

the total amount of traffic passing through the centre o? This quantity scales as n 3 , 
following from scaling arguments using basic stochastic geometry. But in fact one can 
compute the constant of asymptotic proportionality (Theorem [5]). 

2. Indeed, does the scaled flow T n /n 3 have a non-degenerate limiting distribution? (The 
answer is yes: see Theorem [7] and Corollary [9j) 

3. Does uniform integrability hold for the sequence of T n /n 3 asn^ oo? If not then there 
might exist a well-behaved limiting distribution but the mean of T n /n 3 might converge 
to a higher value than that of the limit. Were this the case, it could be viewed as a 
kind of stochastic congestion result. (The results of Theorem and Lemma [S] indicate 
why uniform integrability does hold; it is possible to push this further as indicated in 
subsection 13.41 ) 

Section H] provides a comparison by describing analogous results for flows in cities built on 
grids (Manhattan cities). The concluding Section [5] adds some further remarks, and mentions 
possible future research directions. 



2 Making connections in the Poissonian city 

Asymptotic arguments applied to formulae from stochastic geometry indicate the geometry 
of the routes provided by the unsophisticated semi-perimeter algorithm described above, in- 
cluding their random variation in length, and ways in which they differ from true network 
geodesies between source and destination nodes. We begin by discussing the asymptotic dis- 
tribution of the location and extent of the maximum lateral displacement of a semi-perimeter 
route from the corresponding Euclidean geodesic. 

2.1 Maximum lateral displacement 

Consider the height and location of the maximum lateral displacement of one of the dC(p~ , p + ) 
semi-perimeter routes from a source p~ to a destination p + . Figure H] illustrates 1000 real- 
izations of such routes, with maxima marked by disks, when source and destination are 
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Figure 4: A plot of 1000 semi-perimeters of cells dC(p~ ,p + ) based on a distance n = 
dist(p~ ,p + ) = 1000. The dots indicate the maximum lateral displacements from the horizontal 
axis between source and destination. The figure has been subjected to vertical exaggeration by 
a factor of y/n/A. 



separated by distance n = 1000. Such simulations suggest a limiting distribution under scal- 
ing for the extent and location of the maximum lateral displacement, and stochastic geometry 
arguments show that this is indeed the case. 

Theorem 1. Consider two points p~ = o = (0,0) and p + = (n, 0) located along the x- 
axis, and also a path between these points based on dC(p~,p + ) fl {(x, y) : y > 0}. Locate 
the maximum lateral displacement of dC(p~ ,p + ) fl {(x, y) : y > 0} from the x-axis (and 
thence from the Euclidean geodesic between p~ and p + ) at (nU n ,y/nV„). Then the joint 
distribution of (U n , V n ) has the following weak limit (U, V) as n —> oo; the scaled location U 
is uniformly distributed over [0, 1] and conditional on U = u G (0, 1) the scaled displacement 
V is distributed as the length of a A- dimensional Gaussian vector with variance 2u(l — u). 

Proof. Let p~ , p + be located at o and (n, 0) on the x-axis and let the maximum displacement 
be located at (nU n , \JnV n ) as in the statement of the theorem, so that 

\pnV n = max{|/ : (x, y) G dC(p~ , p + )} 
(nU n ,VnV n ) G dC(p~,p + ). 

(Almost-sure uniqueness of U n is a consequence of the fact that a stationary isotropic Poisson 
line process almost surely contains no horizontal lin es.) 



The proof is a variation on ideas in the proof of lAldous and Kendall! (120081 . Theorem 3). 
Consider the point process formed by intersections of lines l~ , £ + from II subject to the 
following additional requirements: 

1. No further lines from II separate the intersection t~ fl £ + from the segment of length n 
formed between the pair of points p~ = o, p + = (n, 0); 

2. One of the intersecting lines i~ has positive slope, the other £ + has negative slope, and 
neither line intersects the segment formed between the pair of points p~ , p + . 



Topological arguments show that there must be just two points in this point process, one 
above and one below the x-axis, and the point above the x-axis must be located at (nU n , \/nV n ) 
The intensity of the point process in the upper half-plane is given by 

p(x,y) = 

- (sin a + sin/5 — sin(a + /?)) x exp (v ' x2 + ^ 2 + V ( n ~ x ) 2 + V 2 ~ (2) 

where a, (3 G (0, it) are the interior angles at o and (n, 0) of the triangle formed by (x, y) and 
these two points. Here the exponential factor is contributed by requirement [1] above, since 
Slivynak's theorem shows that the unit intensity stationary isotropic Poisson property is 
preserved by conditioning on two lines from II intersecting at (x, y) and then removing those 
two lines. Employing the fact that the intensity of the point process formed by intersections 
of lines from the unit intensity line process II is ~, requirement [2] can be shown to lead to 
the factor 

7T 1 f a [ p 1 1 

— x — / / — sin(0 + ip) dtp d9 = — (sin a + sin/? — sin(a + p)) . 
2 T Jo Jo 2 4 

It follows that (U n , V n ) has joint density on the upper half-plane given asymptotically for 
large n when < u < 1 and v > by 

n 3 / 2 p(nu, y/nv) = 



n I / „ v 2 , v 2 
^^(sina + sin/3 — sin(a + (3)) x exp I — — I y u 2 H h y (1 — n) 2 H 1 

n 3/2 / 1 t> 2 x 



exp I — -— j- I x (sin(a)(l — cos(/3)) + sin(/3)(l — cos(a)) . 



4 V 4m(1-m) / 
Converting the sines and cosines to expressions in u and v , as n — > oo so 

n 3/2 



-p(nu, y/nv 




This can be identified as the joint density corresponding to the limiting distribution of £7 n 
and V n given in the theorem; weak convergence follows from Fatou's lemma. □ 

Simulation studies (from which Figure |4] was derived) confirm these asymptotics. 
2.2 Random variation via growth processes 

The direct stochastic geometry method is highly effective for computing detailed asymptotics 
of mean-value quantities, but leads to burdensome calculations for second-order quantities. 
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Therefore we turn to an alternate approach based on random growth processes. The max- 
imum analyzed in Section 12.11 occurs at the point of intersection of the trajectories of two 
independent growth processes, together representing the lateral deviation of the path from 
the Euclidean geodesic between source and destination nodes. One growth process is viewed 
as starting from the source node p~ and one from the destination node p + , tracing out the 
relevant semi-perimeter of C(p~ ,p + ) by describing the height as a function of arc-length along 
the semi-perimeter. The two processes {Hf : s > 0} are given by heights Hf above the y- 
axis at arc- length distance s along the respective path from the originating node (p~ for H + , 
p + for H~). Let be the angle made by the path with the x-axis, where the angle mea- 
surement is oriented depending on the label "±" so that Gq = n (since the path commences 
by setting out in the opposite direction to that of its goal): then -^Hf = sin 6^ (except for 
isolated points at which the slope of dC(p~,p + ) changes). Changes in occur when the 
path is intercepted by a line from II which also intersects that part of the x-axis with p T 
deleted which does not contain p ± . Indeed, starting at arc- length s along the path from p T , 
the angle Q^ 1 remains constant for an Exponentially distributed length of rate |(1 — cos O^), 
after which the angle jumps to a lower value Qf_ + AB^ < O^L, where the hitting properties 
of Poisson line processes can be used to show 

1 COS (J) 

P \-ASf <(j)\ef_ = e] = £ where < < . (3) 

L J 1 - cos 6 

(Jump processes such as are taken to be cddldg, so that we have lim s ^O s = t , and 
we write 0_ for the process of left-limits.) We suppose the growth processes H evolve for 
all time according to the dynamic described above. Let Xf be the distance from p T when 
resolved along the axis from source p T to destination p ± , so 4-Xf = secO^ (except for 
isolated points of discontinuity of 9). Figure illustrates the H + construction. 



o 

start p~ = o goal p* = (0,n) 



Figure 5: Illustration of H + construction. The growth process Hf tracks height as function 
of arc-length s. The angle of slope G^" is an auxiliary process governed by a Poisson stochastic 
differential equation, jumping when the path is intercepted by a line from II which also intersects 
the negative x-axis. 



Analysis of growth process using Levy processes 

For the sake of clarity of exposition we now drop the superfix ±. 

It is convenient to apply a random time change t = s — X s using the excess of arc-length 
over distance travelled towards goal; in the new time scale the angle changes according to 
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a Poisson process of incidents of rate | while ^ = sin 6/(1 — cos 6), = cos6/(l — cosG). 
This gives a stochastic differential equation for H and X, driven indirectly by a half- unit-rate 
Poisson counting process N via the auxiliary process 0: 

, rr sin 6 

dH = - -dt, 

1 — cos B 

cos 6 

dX = -dt, 

1 — cos B 

dB = ABdiV, (4) 



The distribution of the jump AB is given by (j3J), and the jump AB at a jump of N is 
conditionally independent of the past given B_ at that time. Note that (jlj) can be viewed as 
driven by a marked Poisson process which is obtained by marking the incidents of N by the 
jumps of B, with mark distribution (conditional on the left-limit ©_) given by (J3]). Using 
this terminology and approach, we will now show that the excess <j(n) = inf{t : X t > n} = 
infjt : J* il°ofe ( ~ Lu — n } a ^ Si ven distance X = n has standard deviation asymptotically 
proportional to ^/\ogn for large n. 




Ti 



Figure 6: Illustration of initial segment of a path used to provide an upper bound on the 
excess acquired before X = 0. This initial segment is made up of three line segments, of lengths 
T\, T2, and finally a length bounded above by Xisect/. Here Ti, T2, U are independent with 
distributions given in the text. 



To establish this result it is simplest to consider the growth process begun with Xo = 
and ©0 lying in the range (0, 7r/2]. We must therefore control the amount of excess required 
to achieve this. Since X = at both ends of this segment, this can be bounded above by the 
length of the initial segment of the path indicated in Figure [6j This path uses the following 
directions till it first hits the j/-axis: it first runs along the negative x-axis till it encounters a 
line of angle between | and |; then moves upwards along this line till it encounters a line of 
angle between and |, then along that line. The first segment is of length Ti, Exponentially 
distributed of rate |. The second segment is of length stochastically bounded above by 
T2, Exponential of rate |. The final segment is of length stochastically bounded above by 
T\ sec U, where U has density -7= cos u for < u < | (this uses a stochastic monotonicity 
argument applied to the conditional distribution of U given the angle of the second line 
segment). We may take Ti, T 2 , U to be independent. 
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The quantity Ti + T 2 + Ti sec U is an upper bound on the path length of the initial segment 
and has finite mean given by 8(1 + and finite second moment given by 32(3 + -^(tt + 
log (2 + \/3))). It follows that the contribution of the actual initial segment to the mean and 
variance of the excess is bounded and may be ignored if we can establish logarithmic increase 
in variance of the remainder. Accordingly we may suppose our growth process begins with 
X = and with O lying between and |, distributed according to the density cos^ for 
< 9 < |. (This is because the growth process will intersect the positive y-axis in the first 
intercept which makes an angle with the horizontal in the range (0, |).) 

We now address the question of the variance of u{n) (based on X = 0) in three stages. 
First of all we use trigonometry and probabilistic coupling to relate the negative log angle 
— log G to a Levy subordinator £. We then state and prove a lemma on analogous mean and 
variance asymptotics for r(n) = inf{£ : J* exp(2£ u ) du > n}. Finally we state and prove a 
theorem which uses approximations and coupling to establish the required mean and variance 
asymptotics for a(n). 

For the first stage, note if < < | then by trigonometry and calculus 



2 5 cos 6 2 

62 ~ 6 ~ 1 - cos 6 ~ 6 



< — — < — ■ (5) 



Applying (j3J) to a jump AO = — ©_ produces a unit Exponential random variable 

, (1- cos(-A©)\ 

J = -log — K ——L . 6 

\ 1 — cos B_ / 

Thus we can independently mark each jump of the Poisson process N using the unit Expo- 
nential mark distribution of J . The mark J can be used together with B_ to construct the 
actual jump of O, using ([6]). Bearing in mind that A0 < 0, we may write 

j = /(io g e_)-/(io g (-Ae)) 

where f(x) = log(l — cose z ). Calculus shows that /' > and /" < over the range 
(— oo, log |). Hence if x < log6_ < log | then 

77 < /'(lo g e_) < f'(x) = f smx < 2 = lim 



2 ' 1 — cose^ x-^-oo 1 — cose^ 

But < -AG < 6_ and J = /w^e) f'( x ) dx > therefore 

/'(loge_)0oge_-log(-A6)) < J < 2(loge_-log(-A6)) . 
Hence, for < -A0 < 0_ < f , 

-log(l-exp(-fj)) < -log (l-exp (~ f , {l f gQ) )) < 

< -log^ e - + A9 ^ = -Alog9 < -]og(l-exp(-lJ)) . (7) 

Note for future reference that — log(l — exp (— is distributed as the maximum T' V T" 
of two independent unit mean exponential random variables T' and T", and in particular it 
has mean | and variance |. On the other hand — log(l — exp f— -JTJ) has probability density 
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|(1 — e x )i l e x for x > 0,and the square of its negative exponential (l — exp (—-jX) 
(which will play a role later on) has probability density 

-(l-v^) 1-1 ^ forO<x<l. (8) 

Thus a coupling construction indicated by the inequalities of (GO) permits approximation 
of the negative logarithm of the angle process by rj < — log(0/Go) < £. Here rj, £ are non- 
decreasing pure-jump Levy processes (hence subordinators) which have jumps at the same 
times as those of — log G (namely at incidents of the Poisson counting process N of intensity 
|) but with jump distributions given respectively by the distributions of — log(l— exp (— -JTJ) 
and — log(l — exp(— J /2)). Also for future reference, note that the Laplace exponent = 
-f logE [e~ 9?t ] of £ can be computed as $(<?) = 2(1 ^2+ g ) for ? > _1 > while M * = 6 - f s 
defines a martingale, as does M 2 — \s. In particular M is an L 2 -martingale. 

The discrepancy between coupled jumps of £ and — log can be controlled by 

Now if < a < | then we can use cos 2 | > |, sin | < |, and convexity of tan 2 | over this 
range to establish the simple bound 

1 1 — cos a 1 — cos a 1 tan 2 f la 2 

= < ^ = - H < - H 

/'(log a) a sin a sin a 2 2 2 4 

and so (using 9 < | to apply the above inequality) 



log 



l -exp(- t 7/f(log9_)) 
1 - exp (-\J) 




< l()g|l + 7T~\ "I < f 1 - — < — < -<>xp(-2//_) . 

1 expd^-l 2 J ~ exp(ij) ' 1 1 1 ' 

Applying this upper bound on the jumps, it follows that we can control the total discrep- 
ancy between £ and — log by 

< &-(-k>g(e t /e )) < i^exp(-2^_). (10) 



2 

w<t: 

AA r m >0 

However the right-hand side is increasing in t and has a limit which can be expressed in terms 
of a simple perpetuity. Indeed 

(oo k 
i + Yi n ex p (- 2A -^) 
k=l m=l 

where A m r] is the m th jump of rj. First and second moments of the perpetuity final value 
Uoo are easily expressed in terms of first and second moments of exp (— 2A m r]) (see for 
example IVervaatl . Il979l . Theorem 5.1). However we will require control of exponential 
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mome nts E [exp(zL/ CJO )] for positive z. lAlsmeyer. Iksanov. and Rceslerl (120091 ) and iKellerer 
(119921 ) give results for the general case; however in our particular case the perpetuity mul- 
tiplier exp (— 2A m rj) is positive and bounded above by 1, moreover its probability density 
(El) is bounded above ne ar 1, and thus it follows from monotonicity and the methods of 
Goldie and Griibell fll996l . Theorem 3.1) that E [exp (2: [Too)] < 00 for all positive z. (In fact 
the upper bound re quirement on the probability dens ity can be replaced by comparability 
with a Beta density: iHitczenko and Wesolowskil . l2010l . §4.) 

We can now improve ([5]) to provide bounds in terms of Levy subordinators: 



-2 exp (2f - U a 



5^2 5 cos9 
6 ~ 02 ~ 6 ~ 1 - cos 6 
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Accordingly we first establish a lemma which works with j exp (2£ s ) ds rather than X: 
Lemma 2. Define r(n) in terms of J exp(2£ s )ds by 

r-r(n) 



n 



r-ryn) 

I exp (2£ s ) d s . 
Jo 



Then 



r(n) 



and as n — > 00 so 



log n - 2M r( „) + log ( exp 



|logn + 0(l) 



2 ^r(n) 



n 



3 

§logn + 0( v / log 



E [r(n)] = 
Var [r(n)] 

Proof. First note the trivial estimate establishing finiteness of E [r(n)]: 



(12) 



(13) 
(14) 



n 



exp(2^ s )ds > r(n) . 



The representation (fl2l) was motivated by heuristic time-reversal arguments but is essen- 
tially tautologous: write 



exp (2£ T(n) ) = exp (2M r(n) + §r(n)) 



n x 



exp (2£r(n)) 



take logs, and re-express in terms of r(n). Taking expectations, we then obtain 



E [r{n)\ 



log n + E 



log 



exp (2£ T („)) 



n 



the expectation E [M r ( n )] vanishes because r(n) has finite expectation and so the stopped 



martingale M sAr ( n ) is L 2 -bounded by E 
Consider the variance of rin 



<|E[r(n)]. 



2 I ex P( 2 ^r(w)) \ . ugm g £ 2 _ mar tingale theory we find 



Var 



' 2 /exp(2^ r(n) ) 

r(n) - - log 

3 V n 



Var[M r(n) ] = fx §E [r(n)] 
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So a uniform bound on the second moment of log ( exp ( 2 ^ T ( n > J w [\\ complete the proof. 



To this end, note that Z n = exp(2£ r(n )) constitutes a Lamperti transformation of the 
subordinator 2£ and therefore defines a self-similar Markov process Z. Using the scaling 
property for Z, 



E 



exp(2£ r(n) )\ 2 exp(2£ r{n) , 
log I ; — < 1 



n 



n 



E[(logZ 1 ) 2 ;Z 1 <l|Z = i] 

< 4E[Z^\Z = ±] 



and iBertoin and Yorl (120051 . Formula (20)), drawing on iBertoin and Yorl (]200ll ). can be ap- 
plied together with the calculation of the Laplace exponent <3>(g) of £ to show that 



(1 + {n - l)e~ n/2 ) < |(1 +2e 



For later use we note that each integral moment of Z 1 1 is bounded, since for p > 

2 P I / r nP / 2 



2p+l 



n P e ~nP/2 



)P -l e -v/2 d , 



(15) 



To bound E 



log 



cxp(2g T(n) ) 



cxp(2g T („)) 



> 1 



r' n = inf{t : 2£ t > log ft}. Hence 



for n > 0, note that r(n) < tL + 1, where 



E 



log 



exp(2e 



r(n). 



exp(2£ r(n) ) 



> 1 



< E 



foe i \ 2 exp(2£ r(n) ) ' 
(2£r;+i - logn) ; > 1 



< E 



- log ft)' 



Now 2£ 7 j + i — log ft is the independent sum of a summand of distribution 2£i and a summand 
which is the overshoot 2£ T ^ — log ft. Since the jump distribution of 2£ is 2(X" V T") (for two 
independent unit exponential random variables T' and T"), it follows by the memoryless 
property of the exponential distribution that the overshoot is some mixture of the distribu- 
tions of 2T" and 2(T" V T"), hence stochastically dominated by 2(T' V T"). Consequently it 



follows that E 



(2£<+i - log ft)' 



< 315 
— 16 ' 



Thus E 



log 



exp(2g r(rt) ) 



< |(1 +2e 2 ) + -j^ ; and hence the lemma is proved. 



□ 



Note that upper bounds for E 



log 



exp(2£ T(n) ) 



can also be obtained using the techniques 



of IBertoin and Yorl (120021 ). but these bounds diverge to infinity with n. 

We are now able to attack the asymptotic behaviour of mean and variance of a(n): 

Theorem 3. With a{n) defined as above, so that n = X a ( n ) for n > 0, for large n 

E[a(n)} = f logn + 0(1), (16) 



Var [cr(ft)] 



§ logn + 0(1), 
§logft + 0( v / log 



ft 



(17) 
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Proof. Lemma [2] provides partial control on the asymptotic mean and variance of cr(n) via 
ATI) , since X t < J* exp(2£ s ) d s and so r(-^-n) < a(n). By the representation (Tl2"j) and 
integration (since 9 has density cos 6* over (0, |)) we find that 



E 



'r(fn) 



f logn + ~E 



0(1) 



Var 



r(fn) 



log 



exp(2^ r(0 



9gn 



20 
27 



logn + 0(l), 
logn + 0(1) . 



The second moment of log 



/cxp(2g T(e 2 n) ) 



20 
27 



being bounded, it follows that Var 



r(fn) 



logn + 0(\/logn). 

So consider the lower bounds from fflTT) . delivering an upper bound for o~(n) via: 



02 
u o 



6 



_exp(2£ s -[/ s )-~ ds < / — - - ds < X, 



©2 



6 



First observe what happens after the stopping time given by 



1 + — eg) x inH f : 



e? 



> 



2 
u o 



+ 



We find that 



e 



+K 



> 



2 el 
92 92 

w w s+ft 



5_ ©?. 



6 



5 

> 

6 ~ 



2 el ^ 

^0 ®S + K 



and moreover 



JH lds 



> 



1 + lC)2 
1 + 12 U 



Aft2 
12^0 

1 + 12 U 



2k 

02 
U 



> 



This will allow us to disregard the effects of the — | term, so long as we can bound the second 
moment of k. To this end introduce (, a Levy subordinator jumping at the incidents of the 
underlying Poisson process N, with jumps coupled to those of the other processes so that 
C < V < — log(©/©o)j and with 



AC = Ir 2 

1 -log(l-exp(-- 1 7))>l 



Since these jumps are all of sizes or 1, we can view ( as a Poisson counting process run at 
rate v < \. Moreover, since < O < we have 



2 • 



K 



< 



5 

12 



0^ inf U : 



Q I exp(2C 4 )> Q2 



2 5 
+ 6 







< 



1 + 



5^ 
48 



inf U : (t > 7T log 1 + 



57T 2 

48 



and the boundedness of the second moment of the right-hand side of these inequalities follows 
directly by comparison with a Gamma random variable, derived from basic Poisson process 
properties. Hence E [k 2 ] < 00. 
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So consider £ s = £ K+S — £ K , and r(n) = inf{t : f* exp(2£ s ) ds = n). We will bound ir(n) 



above by a stopping time K + r(^-n) + p, where p is chosen to compensate for the undershoot 
of n at time t = k + r(-^-n) caused by the U contribution in 



/ (<r xp(2f - 



u. 



5 Id, 



< X f . 



We find 

"K,+r{-^-n)+p 



2 

02 



B 2 



ds > 



r(—n)+p 2 







2 exp(2£ s - 



> 



exp 



-U 



b 2 



n 



Thus we can choose p to compensate for the undershoot so that 



)ds 

J 2 -pexp(2^ 



exp 



-E7 



02 



^+ 0|Pexp ( 2£ 



e 2 



this is fulfilled by the choice 



P 



exp I U 



e 2 

K+r(-^-n) 



- 1 



/exp(2£ 



n : 



-i 



x 



e 2 / 



(■-)'-■ 



V 



/ 



Now the first factor is bounded above by exp(C/ 0O ) , and we have already noted that E [exp(^?7 00 )] < 
oo for all z > as a consequence of perpetuity theory. The second factor is distributionally 
a randomization over m of (exp(2£ T ( m ))/m) _1 , and we have already noted that Lamperti 
transformation and the results of Bertoin and Yor allow us to bound E [(exp(2£ T(m ))/m)~ :P ] 
uniformly in m for any fixed p > 1. Thus we may deduce E [p 2 ] < oo as a consequence of the 
Cauchy-Schwartz inequality. 
Accordingly we find 



e 2 



a(n) < k + r(-^-n) + p 



bounds cr(n) above by a stopping time with expectation |logn + 0(1); moreover we may 
apply the representation ffT2l to deduce that up to terms of 0(1) second moment 

« 2M on -2M 



r(f n) - r(f n) 



B 2 

~/ U . 



k+t(— n) 

which itself must be of uniformly bounded second moment: 

2" 



e 2 



E 



2M 



e 2 



-2M 



K+r(-^-n) 



e, 



E 



K + T' 



e 2 



, C)2 



: ((E [«] + logn) - (logn)) + 0(1) 



0(1). 



Since 



'6 2 



< 



~ © 2 

cr(n) < K + r(^n)+p. 
16 



~ 2 2 

and and differ only by a quantity which has 0(1) second moment: 



E [<r(n)] 
Var [cr(ra)] 



E 
Var 



r(%n) +0(1) 



O A /Var 



r(%) +0(1). 



Consequently the theorem is proved as a consequence of the upper bound asymptotics ( fT8l) 
established at the beginning of this proof. 

□ 



A Brownian digression 

Because we can compute Var \ M t ] = |t, we can use martingale central limit theorem ideas 
flRebolledol . Il980l : IWhittl . 120071 ) to show 

t « | ^logX T + v 7 ^ + 20 - log2 - log^ exp ^-§u + a/P^ dw^) , 

for i? a standard Brownian motion not independent of B, and a T a stopping time for B with 
expectation E [a T ] = t. The distribution of the Dufresne integral 



exp -iu+ \ \B T U du 



is known explicitly (jDufresnd . Il99dl ; lYorl . Il992l ); however its contribution to the above is 
dominated by other terms. 



Recovery of logarithmic excess result 



Inspection of the growth process analysis shows that the logarithmic excess occurs before 
(say) time n/logn, whereas our discussion of the maximum lateral deviation of dC(p~,p + ) 
(with dist(p~ ,p + ) = n) shows that the intersection of the two growth processes occurs at 
x-coordinate uniformly distributed over the range from p~ to p + . 

Hence the asymptotic excess of the upper or the lower semi-perimeter route for dC(p~ ,p + ) 
must have l eading term | log n, in con cordance with the asymptotic for the mean excess 
obtained in 
obtained by 



— — 1 J 1 

Ald ous and Kendall! (120081. Th eorem 3) and similar higher- dimensional results 



Renyi and Sulankd . Il968l . Satz 5). 



Boroczkv and Schneider! ( 2008 . Theorem 1.3) (compare the planar arguments of 



2.3 True geodesies 

We can now deduce that the two semi-perimeter routes provided by dC(p~ ,p + ) will often not 
be geodesies; the boundary dC(p~,p + ) is composed of the initial parts of four independent 
growth processes, contributing four independent initial excesses each of mean |logn and 
variance proportional to logn; the remainder of the excess will be of order less than y/\og n. 
Accordingly there is an even chance that the least excess is achieved by crossing over from top 
side to bottom side so as to use the smallest possible | logn ± const, x ^f\ogn contribution. 
Calculations of the caricature of Section 11.11 makes it plain that such a crossover can be 
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achieved at the very modest price of adding just a bounded term to the excess, and therefore 
there is a substantial positive chance that one of the crossover routes is shorter than either 
of the semi-perimeter routes. 

In fact we conjecture that the two semi-perimeter routes provided by dC(p~ ,p + ) are never 
geodesies; in particular it should be possible to achieve modest reductions in the excess by 
using crossovers very close to source and destination nodes p~ and p + . 



Lower bound for the Poissonian city 

Nevertheless the semi-perimeter routes supplied by dC(p~,p + ) are good approximations to 
true geodesies; we show this by establishing that their mean excess can be compared with a 
lower bound on possible path-lengths. Indeed, because we are working in the specific situation 
of a Poisson line pr ocess we can derive a strong er and simpler version of the Q(y/logn) lower 



bound argument of lAldous and Kendalll (120081 . Theorem 2). 



Theorem 4. In the Poissonian city network, consider any path from p to p + (in the sense 
described in Section ing) . If dist(p~ , p + ) = n then the path must have mean excess exceeding 

log 4 — - J log n + o(log n) = 0. 136294 . . . log n + o(log n) . 

Proof. Let C(o, +) be the cell containing the positive x-axis of the tessellation formed from 
the Poisson line process IT by deleting all lines intercepting the positive x-axis. Consider the 
vertical line £ x through (x, 0), and let — L~, be the distances along this line to dC(o, +) 
running down and up. Any network geodesic 7 from o to any other point p on the positive 
x-axis, constructed according to the recipe in Section 11.21 must lie between locations — L~ 
and L + on £ x . This is a consequence of the convexity of C(0, +). Consider such a geodesic, 
or indeed a general regular path 7 lying within these bounds, and let 9 X be the angle made 
with the horizontal by 7 when encountering £ x for the first time. If dist(o,p) = n then the 
mean excess of 7 must exceed 



E 



pn 1 rn rn PTr/2 

/ (sec^-l)dx >\j E[^]dx= / / P [\0 X \ > u) udu dx 
Ji J Ji Ji Jo 



Consider the probability of there being no lines of II which both (a) hit £ x in the range 
between — L~ and L + signed distances from x-axis and (b) form an angle to the horizontal 
which is less than u in absolute value. The density of the angle to the horizontal is | cos 6 
for — |- < 9 < It, while the patterns of lines hitting £ x above and below the x-axis are 
independent. Consequently 

¥[\6 x \>u] > E [exp(-(L- +L+)sin M )] > (E [exp {-uL+)]f . 

Considerations from stochastic geometry show that 

E[exp(-wL+)] = f P [exp > z] dz 

Jo 

/"OO /"OO 

= / P [exp (-uL+) > e~ s ] e~ s ds = 1- e~ s P[L+>f]ds 
Jo Jo 





-»-!( 


L exp ( 
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Consequently 



E 



(sec 9 r — 1) dx 



/n fir/1 
J F[\9 X \ > u}udu dx > 

n fir/2 / poo 
> I I II- 



1 JO 







exp | — s — \ [ \l x 2 + 7?z — x \ \ ds\ udu dx . 



Suppose that in the above we could replace y x 2 + ^ — x by its upper bound \s 2 / (xu 2 ). 
We would then need to estimate 



n fir/2 



1 JO 



I — exp — s 



Au 2 



x 



ds\ udu dx . 



Now we can estimate 



Pexp \ J exp ^-y ) ds 



usin g classical re sults on Mill's ratio, namely ISampfordl (119531 ) excellent upper bound (see 
also IBarica . 120081 for a treatment based on a monotone form of l'Hopital's rule): 



p 2X rx 

eX P I 17 

^ / .ip 



s 2 \ 4 

exp I ds < — . 

2/ v /^g + 3 p 



for p > — 1 . 



(19) 



Accordingly, setting v = y/xu, as n tends to 00 so 

»n /-ir/2 

ds ) udu dx 



1 - / exp f -s - 

1 Jo \ Jo v 

yfinr/2 



r-n r-Jxir/2 / r-oo W 

k Jo V~l ^{- S -^) ds J Vdv 

II I ; Vdv 

Vv 2 + 4 + 3v x 



d x 

x 



1 JO 



1 JO 



\/v 2 + 4 — d \ dx 
, vdv — 

Vv 2 + A + 3v x 



(log4-f)log 



n . 



The proof is completed by bounding the error arising from the approximation of -\f x 2 + s 2 /u 2 - 
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x by \s 2 1 [xu 2 ]: 



n I"k/2 / poo 

II—/ exp ( — s — \[ \l x 2 + ^2 — x ) ) d s ) udu dx 

' I •/(! V JO 




an/2 



1 — e s [ exp 



l(\/^ + ^- S ))- eX P (-i^)) ds 



oo 

2 







»n /-Tr/2 



d s ] udu dx 



L L ex p(- s -ii) ds ) wdudx 

x / e~ s (^exp ^^xM^-x^ -exp ^-JL_JJ ds«dH dx. 



We bound the second term by invoking iBirnbauml (jl942j)'s very good lower bound: 



;°° / s 2 \ _ ^ vOT-v^a 

1 — / exp — s — d s < 



o 



— p ds < - . 

4M 2 X J y/xU 2 + 1 + v/xii 



Thus it suffices to bound 



n rir/2 — /-oo / -±IJ X 2 




, v / e e V / — e dsudu di 



1 JO JO 

A JO V 1 +^r+ 1 Jo I J u x 

1 /-/2 r -2 ! / r— \l dx 

< - / / / e mm < 1, — — \ x 2 + % - x ^ ds dw — = 

~ 2Ji Jo Jo I 2 V V J ) 



X 



where we use the fact that over the range < p < oo the function p t— > e~ p is nonnegative, 
decreasing, and Lipschitz with Lipschitz constant 1, and also that 1 + — y/1 + p < |p 2 if 
< p < oo (use finite Taylor series expansion). 

We split the x-integral at x 3 = s 4 /(16m 4 ); the first part is bounded by 



1 r r /2 r (s4 ^ wu ^ 1/3 _ s r i s 4 i dx , , 

- e min < 1, — > — d« ds 

2 ./n ./n ./n 1 16 K 4 X 3 f 



JO JO 



■y roo rn/2 p(s i / (16m 4 )) 1 / 3 j 



2 




Jo Jo 



/ e- s ^Lduds = Ivr^rd), 

Jo V% 
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while the second part is bounded by 



i r° r /2 r r 1 s 4 i dx , , 

e mm < 1, — — > ^ d« ds 



27 7o A^ 4 /(i6« 4 )) 1 /3 " " I' 16m 4 x 3 J 

- /°° /°° e - s -4^^dnd S = ^Tr^rd). 

2 Jo Jo ./(^/(leu*)) 1 ^ leu^V^ 

This establishes the desired result. □ 

The results of this section justify the focus in the remainder of this paper on the semi- 
perimeter routes provided by dC(p~ ,p + ): while semi-perimeter routes do differ from geodesies 
nevertheless their use does not incur a great penalty; they are produced by a geometric 
algorithm which is certainly unsophisticated but on the other hand is explicit; and they are 
amenable to exact calculations. 



3 Traffic flow in the Poissonian city 

We now consider traffic flow in the network produced by this Poisson line process, and to do 
this we first compute the mean flow through a line at the centre of the disk. More formally, 
we condition on there being a (horizontal) line of the line process running through the origin 
o, and consider the flow through o which results if every pair of x and y in ball(o, n) generates 
an infinitesimal flow of amount dx dy divided equally between the two possible routes given 
by the semi-perimeter algorithm (see Figure [7]). 




Figure 7: Illustration of the flow generated between two points in a "Poissonian city" . 
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3.1 First-order calculations at the centre 



Recall Equation ([I]) from the introduction: the flow through the centre is measured by the 
4- volume of T> n where 

£>n = {(p~,P + ) e ball(o,n) 2 : Pi<Pi, o G dC(p~,p + )} , 

and we seek to understand the large-n statistical behaviour of this 4- volume. Indeed (bearing 
in mind that we have conditioned on there being a line through o) the distribution of the 
total traffic through the origin o is given by the distribution of 



T n = \ if l[(p-, P +)ev n ] dp dp + = \ / / \- <p +,oedC(p-, P +)} d P d P" 

J J J ^ball(o,n) 2 

The mean can be obtained asymptotically using direct arguments. 
Theorem 5. The mean flow through a line at the centre is given by 

pir pn pn 

E [T n ] — / / / exp(-i(r + s-p))rdr sds9d9, 
Jo Jo Jo 



(20) 



(21) 



where p = \Jr 2 + s 2 + 2rs cos 8. Asymptotically as n — > oo ; 

E [T n ] ~ 2n 3 . (22) 

Proof. Equation ff2T]) follows from simple stochastic geometry of Poisson line processes (illus- 
trated in Figure [8]), taking care not to double-count flow between the un-ordered points p + 
and p~ . Note that when p + and p~ are on opposing sides of the line conditioned to hit the 
origin then none of the flow between these two points will run through the origin. Indeed, 
mean flows between points p + , p~ in the upper half-plane will account for exactly half the 
total mean flow through the origin. 





^^^^ r 




/ s \ 




9 V 


/ ° 





Figure 8: Illustration of the geometry represented by the multiple integral in (|2ip using r, s, 6. 
The segment p~p + is not separated from the origin o exactly when no lines of the line process 
pass through both of op~ and op + . 



22 



The asymptotics follow by application of scaling by n, the symmetry between s and r, 
and the inequality y/l — 2z < 1 — z for z > ~. Indeed . 



pir pi ps 

E[T n ] = 2n 4 / / exp(-f(r + s-p))rdr sds9d9 
Jo Jo Jo 

= 2n A J exp (-f (r + 1 - Vr 2 + 1 + 2r cos 0)) r dr s 3 ds #d#. 



The region of the integral corresponding to J^ 2 J" J" is easily seen to be bounded above 
by 6tt 2 ti 2 , while the region Jg is bounded above by 87m 2 /(l — cos# n ) when 9 n > 0. 



Consider the region J Q n J Q J" . Using a Taylor series expansion of 1 — y/1 — 2z, and the 
approximation 0/ sin 6 \ 1 as 6 \ (so long as < 9 < tt/2), we deduce that 



EfTJ ~ 2n 4 f f f exp ( , nSr s u] dus 3 dsrdr 

1 1 Jo Jo Jo V 2(1 + r) ,/ 



= 2n 4 Z" 1 / >l 2(1 + r) g 3 d g rdr = 2n 3 . 
Jo 7o nsr 

□ 

Taking some extra care over the analysis, it is possible to bound the order of the error in 
Equation (122]) ; we state this without proof. 

Corollary 6. The asymptotic in Theorem^ can be sharpened to 

E [T n ] ~ 2n 3 + 0(n 2 y/n) as n -> oo . 

3.2 Mean flow averaged over entire disk 

Regional variation of expected flow over the disk is to be expected: flow at the boundaries 
should be lower than at the centre. Indeed, one can calculate the mean flow per unit length 
in the network as follows. 

The mean total length of the intersection of the Poisson line pattern with the disk is given 
by (mean number of lines hitting disk) x average intersection length: 

(2ixn) x — / 2Vn 2 — x 2 dx 

and this is therefore the mean total network length in the disk. 

On the other hand, the mean Euclidean distance between two independent uniformly 
random points in the disk is given by 

— - / / — - / / *J u 2 + v 2 — 2uv cos(a — (3) vdv dfludu da 
™ l J Jo nn 2 Jo Jo 



/ / ^u 2 + 1 - 2wcos(#) udu d9 = — / / s 2 dsd 
Jo Jo 5w Jo Jo 



-n . 



5vr J Jo 5ti Jo Jo 45vr 

where the first step uses various symmetries and re-scaling, and the second step changes to 
polar coordinates based at u — 1 and 9 = in an implicit use of Crofton's method. (This 
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calculation is a special case of a classic calculation in geometric probability, surveyed in 



Santalolll976l . II. 12.7 Note (6).) By the previous results on lengths of network geodesies, the 



mean network distance differs only by an extra logarithmic contribution. 

Hence the mean flow per unit length, if each pair of points exchanges just one infinitesimal 
unit of traffic and this is averaged over the network, is asymptotic to 

1 (Tin 2 ) 2 x j^n 128 , , 
-- o 2 45w = n 3 = 1.9052... n 3 . 

2 45vr 

2 

This analysis does not take account of routes which move outside the perimeter of the disk; 
however the effect of these routes can be shown to be negligible. (The key observation is 
based on Theorem [TJ if both source and destination nodes p ± are at least 2^J (1 + e)n \ogn 
from the perimeter of a disk of radius n, and n is large enough, then points outside the 
disk have probability at most 0(n~^ 1+e - ) ) of lying within C{jT ,p + ). Thus mean total length 
outside disk is a boundary rather than an area effect.) In conclusion, and unsurprisingly, 
mean flow over a typical line is slightly smaller than mean flow over a line at the centre of 
the disk. 



3.3 An improper anisotropic limiting line process 

We can represent the scaling limit of the distribution of traffic flow through the centre of the 
Poissonian city by using an improper stationary anisotropic Poisson line process. 

We use the alternate coordinatization of a unit-rate stationary isotropic Poisson line 
process as in Figure El using coordinates x of the intersection along the x-axis and 9 for line 
direction. Then the x-axis intersections form a stationary Poisson point process, while the 
angle density is | sin 9 for < 9 < ix. 




Figure 9: The Poisson line process can be represented in terms of a Poisson process of points 
scattered along the x-axis, through each of which there runs a line making an angle 6 £ (0, it) 
with the x-axis, with density \ sin 6. Here we show this against a backdrop of the disk ball(o, n). 



Re-scale to shrink the x-axis by a factor 1/n, so x = x/n. Guided by previous results, 
shrink the ?/-axis by a different amount 1 / y/n, so y — yj \fn. Thus 9 is transformed into a 
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new angle (see Figure [TU]), where 



nta.n6 



tan0 = 

In the new coordinates of x and 0, the line process can be parametrized as a non- stationary 
Poisson point process on x : <p space with intensity 

.2 JL 1 



sm sec 



2 (l + itan 2 

V n 



.3/2 



dsdi 



sm sec 



We can represent this as a coupling construction: based on an improper stationary anisotropic 
Poisson line process with intensity | sin 0sec 2 dxd <j) in x : coordinates, we can achieve 
a proper stationary isotropic Poisson line process at scale n by randomly thinning the lines 
with retention probability depending monotonically on the line slope. 




Moreover this limiting object may be cleanly represented using a further set of coordinates. 
Represent each line of the line process by its intercepts y + and y_ on the vertical axes x = 1 
and x = — 1 (see Figure [TTI) . Then the intensity becomes 

In particular, while the new improper Poisson line process is anisotropic nevertheless it does 
possess special affine shear symmetries; namely the symmetries produced by those area- 
preserving linear transformations which leave invariant all vertical axes. 
This construction enables us to identify the limiting behaviour for T n : 

Theorem 7. The scaled quantity T n /n 3 has a limiting distribution given by the analogous 
flow at the centre for the limiting improper stationary anisotropic Poisson line process given 
above. 

Indeed, we can relate scaled fmite-n instances to the limiting case by a coupling argument 
involving the addition of further lines; however the resulting almost-sure limit is not monoton- 
ically decreasing since it will involve a double integral (as in (1201) ) taken over ever-increasing 
regions of the vertical strip in Figure [TTJ 
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Figure 11: This illustrates the parametrization of lines from the improper limiting line process 
as intercepts on two parallel y-axes. 



Before proving this theorem we show that the mean flow at the centre for this limit is in 
agreement with the asymptotics given in Theorem 

Lemma 8. The flow at the centre for the limiting improper stationary anisotropic Poisson 
line process given above has mean value 2. 

Proof of Lemma. Consider first the probability that the line segment from (—a, t) to (u, t) 
is not separated from the origin by the improper line process. The mean measure of lines 
implementing such a separation (measured using the intensity measure of the improper pro- 
cess) is A + B + C, where the contribution A arises from separating lines hitting the top left 
shaded triangle in Figure [121 the contribution C arises from those hitting the top right shaded 
triangle in Figure [121 an d B is derived from the contribution of the remaining separating 
lines. 
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Sure 12: Computing the mean flow through the centre for the flow based on the improper 


anisotropic line process limit. 



Then 
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and similarly C=^-(^ — l); finally 



B 



(t + s) ds 



Consequently 



A + B + C 



4 \ a m 



Since we are dealing with a Poisson process, the required probability of the line segment from 
(— a, t) to (u,t) not being separated from the origin is 



, t2 f 1 1 
exp -- - + - 

4 V a u 



(23) 



Consider the special affine shear symmetries which leave the x = axis fixed. Because of 
this symmetry group it follows that the probability of the line segment from (—a, b) to (u, v) 
not being separated from the origin agrees with (1231) when the line segment from (—a, b) to 
(w,t>) passes through (0,t). This occurs when t = b "t atl ; moreover if we set s = b — v then 



a+u 



dt ds = db dv. Accordingly the mean 4- volume of the region representing the flow through 
the centre is given by 



i 

o jo 



exp 



1 r l 
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1 1 
- + - 

a u 



ds dt da du 



o 



'11 
a + u) I - + - 



a 



u 



exp 



-— (- + -Y) tdt da du 
4 \a uj J 

2 / (a + u) da du 
Jo Jo 



-- 2. 
□ 

Proof of Theorem^ Consider the affine shear transformation T n : [—1, 1] x (0, oo) [—1, 1] x 
(0, oo) given by T n (u, v) = (nu, \/nv). Define coupled random functions 

I n : ([-1, 0] x (0, oo)) x ([0, 1] x (0, oo)) - {0, 1} , 

In(P,Q) = I[r„peball(o,n)]I[T„ 9 eball(o,n)]I[oeC(r„p,r n g)] • 

So I n depends implicitly on the underlying Poisson line process: the previously described 
coupling construction shows that we can couple different Poisson line processes for each n so 
as to arrange that I n (p, q) — * I(p, q) almost surely for Lebesgue almost all p, q, where I(p, q) 
is given by an analogous construction based on the limiting improper anisotropic Poisson line 
process, but not using T n . Moreover we can realize T n using 



T n /n 3 = \J J I n (p,q)dpdq. 



From Theorem O and Lemma [S] it follows that 



E 



I J J In(p,q)dpdq -> E \J J I(p,q)dpdq 
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On the other hand if we restrict consideration to the finite measure space fix ([—1, 0] x (0, K))x 
([0, 1] x (0, K)), for any fixed K, then we may deduce /^-convergence of I n to / via the dom- 
inated convergence theorem, since the indicator functions /„ are bounded. (Here (Q, J, P) is 
the underlying probability space.) 

It then follows from non-negativity of the I n , I that we can apply Fatou's lemma to 
deliver /^-convergence on all of Q x ([—1,0] x (0, oo)) x ([0, 1] x (0, oo)), and so can deduce 
convergence in distribution as required: 

T„/n 3 = \J J I n (p,q)dpdq ^ \ J J ' I(p,q) dp dq 



viewed as random variables (functions of uj e Vl). 



□ 



Note that this proof also establishes uniform integrability of the sequence of random 
variables {T n /n 3 : n > 1}. 

It is apparent from this construction that the limiting distribution is largely insensitive 
to modest variations in the geometry of the city (disk ball(o, n), or square of side 2n, or . . . ); 
however we will not explore this here. 

In principle it is possible that the limiting distribution of T n /n 3 might be degenerate. 
That this is not the case follows rapidly from representation of the limit in terms of the 
improper anisotropic Poisson line process. 

Corollary 9. The limiting distribution ofT n /n 3 is non- degenerate. 

Proof. Let Ek be the event 

E k = [ there is a line connecting {-1} x [0, |] to {+1} x [0, |] ] . 

Then Ek has positive probability for the improper anisotropic Poisson line process; moreover 
E 1 , E 2 , . . . form a monotonically decreasing sequence of events whose intersection is a null- 
set. It follows from elementary measure theory that 



l -E 
2 



J J I(P,q)dp dq ; E k 



0. 



However simple constructions show positivity of the conditional probability 



E 



J J I{P^q)dp dq 



Ek 



> 



for each k. Because each event E k is of positive probability, it follows that for each e > we 
can find k such that 



< - 2 E 



J J I(p,q)dp dq ; E k 



< e, 



and this establishes non- degeneracy of the limiting distribution of T n / n 3 



□ 
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3.4 Higher-order moments 

We have established that one can produce a coupling construction to show that T n /n con- 
verges almost surely and indeed in mean value to the corresponding quantity for the improper 
stationary anisotropic line process. In fact it is possible to establish convergence of moments 
of order 2 — e for e £ (0, 2); computations show that the second moment E [T%] is bounded 
by const, x n 6 hence a uniform integrability argument may be applied. The tiresomely com- 
plicated computations are omitted; we simply indicate the general approach. The second 
moment can be expressed using an eight-fold integral 



n pn P7T fti rn rn pit 



JO JO JO JO JO JO JO 



F[E(r,s,9,a,u,v,<f>,P)] 

d(3 d(f> udu v dv dadOrdrsds. 



where E(r, s, 9, a, u, v, 0, (3) is the event that neither of two line segments ((r, a)-(s, a + n — 9) 
and (u, (3, )-(v,/3 + -n — 0) when written in polar coordinates) are separated from the origin 
by the line process (see Figure fT3l). 

Analysis of this eight-fold integral is complicated because the probability 

P[£(r,s,0, 0,0)] 

takes on around eight different analytical forms according to the relative geometry of the line 
segments. Case-by-case mathematical analysis shows that the multiple integral is bounded 
by const, x n 6 . 




Figure 13: Illustration of the event E(r, s,9,a,u,v, 4>, f3), which happens if neither of the two 
red segments are separated from the origin by the line process. 



4 Comparison with a Manhattan city 

It is natural to ask how the Poissonian city might compare with an alternate Manhattan city 
based on a more conventional Cartesian grid structure of roads (a "grid city"). Consider the 
case of a disc of radius n furnished with roads arranged in a fixed unit-length grid structure. 
Suppose we wish to connect from the point at — (u, v) to the point at (x, y) (using Cartesian 
coordinates). If u, v, x, and y are all non- negative then there are a wide variety of possible 
geodesic connections. Working with the most direct analogy to the results described in 
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Sectional suppose that traffic from —(u, v) to (x, y) divides equally between the two extreme 
geodesies running from —(u,v) to (x,y). Working to 0{n 3 ) (allowing us to ignore some 
double-counting), the total traffic through o is made up of four contributions, arising from 
(i) u — 0, v > 0, (ii) u > 0, v — 0, (iii) x — 0, y > 0, (iv) x > 0, y — 0. Each case individually 
contributes a term of the form 



(x,y)>0 
x 2 +y 2 <n 2 

where we omit the negligible contributions arising when there is just one geodesic between 
source and destination. We can sum these contributions, since the effect of double-counting 
is again negligible. Thus the total flow through o is 7m 3 + 0(n 2 ). If we fix attention on total 
flow through one of the bonds attached to o (so as to establish comparability with the results 
of Section [3TT1 for the Poissonian city) then we obtain ^n 3 + 0(n 2 ), compared with mean flow 
for the Poissonian city of 2n 3 . 

However account needs to be taken of the greater total length of the grid network. Mean 

2 2 

total network length produced by the unit intensity Poisson line process is ZL y-, compared 
with 2 x nn 2 for the unit grid structure. Thus a comparable grid structure has to be based 
on segments of length - rather than 1. The flow produced through a bond attached to o by 
such a grid using the above protocol will be of order 



(nn 2 ) 2 



\2\2 



X 



(fn) 3 = 2n 3 



where the second term computes the flow through a centre bond when re-scaling from n to |n, 
and the first factor is a correction to ensure that total traffic is (tcti 2 ) 2 not (Tt(\n) 2 ) 2 . Thus 
the Poissonian city is comparable to this grid structure in terms of mean flow at the centre. 
However the grid structure with this protocol can be shown to have the undesirable feature 
that asymptotically a definite proportion of the total traffic (around 2%) occurs outside the 
disk. 

In fact one can also carry out calculations for a somewhat more demanding situation in 
which asymptotically the traffic stays within the disk as in the case of the Poissonian city; 
instead of supposing the traffic is divided equally between the two extreme geodesies, we 
suppose that it is divided equally amongst all possible geodesic connections. In effect, the 
actual geodesic is chosen uniformly at random, and in that case the probability that the 
resulting geodesic passes through the origin (0, 0) is 

P [Bin( M + y, p) = u] P [Bin(x + y,p) = x] (T) (?) m) 

P [Bm{u + v + x + y,p) = u + x] ' C^ST*) ' 

Here p = ^ but in fact the same result is obtained for any < p < 1. Choosing p = , u+x , 

r 2 ' ° c u+v+x+y 

so as to control the denominator, and applying Stirling's formula, we can use Taylor expansion 
to approximate this by 

1 (u+v+x+yf/ 2 f (u+v+x+y)(uy-xv) 2 \ /ggN 

^ {u+v)(x+y)(u+x)(v+y) V \ 2(u+v)(x+y)(u+x)(v+y) J " V > 

Using polar coordinates based on an axis at 45° to the Cartesian axes, and a Gaussian 
approximation based on sin 2 (# — (f>) w (9 — (p) 2 , we obtain a heuristic approximation for large 
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n for the flow between two opposing quadrants of the disc: 



(u+v\ (x+y\ 

EEE^yy - 2« s - (26) 

u,v>0 x,y>0 \ u+x ) 
u 2 +v 2 <n 2 x 2 +y 2 <n 2 

Conversion of this approach into a rigorous asymptotic ar gument wou l d require close atten- 



tion to detailed asymptotics of the Binomial distribution (jLittlewoodl . Il969l : iMcKayl . Il989h . 
However there is an alternate argument which is more easily made rigorous: the expression 
(1261) can be re-expressed in terms of a symmetric random walk X as 

^2 F i Xu + v =v-u \ X u+v+x+y =v - u + y -x] . 

u,v>0 x,y>0 
u 2 +v 2 <n 2 x 2 +y 2 <n 2 

Under the Z-action u—HL — l,v—*v + l,x—*x + l,y—*y + l, a statistical pivot argument 
applied to the summand generates a probability distribution on even integers or odd integers 
according to the parity of u+v, x + y. An argument using the Hoeffding inequality quantifies 
how this probability distribution concentrates around its mode; the denominator is controlled 
by choosing the probability P [X n+ i = X n — 1} = p = ■ Thus it can be shown that the 

asymptotic behaviour of the quadruple sum is given by the number of Z-orbits containing 
modal representatives close to the line between —(u,v) and (x,y). This number can be 
expressed as a sum susceptible to elementary asymptotic analysis, finally yielding a rigorous 
argument for the asymptotic given in (1261) . 

Consider now the total flow through a unit-length bond I connected to the origin. This 
will equal half the total flow through the origin, which itself can be viewed asymptotically as 
the sum of two equal components from two different pairs of opposing quadrants. Thus this 
total flow is again asymptotic to 2n 3 . 

Again a grid structure comparable to the Poissonian city must be based on segments of 
length - rather than 1, and a scaling argument then shows that such a grid structure produces 
mean flow at the centre which is asymptotic to - x 2n 3 = 2.54648 . . . n 3 . Thus traffic through 
the centre under this protocol is about 25% higher in a comparable Manhattan city. 

Geodesies in the Manhattan city are on average longer than those in the Poisson line 
process; we can in fact argue in a manner analogous to that of section 13.21 to show that 
mean network distance between two independent uniformly random points in the disk will 
be asymptotic to 

nx — / (sin 9 \ + cos 9 ) d = nx 



457T 27T Jq 457T 71 

so the mean network flow over the whole disk for the grid will again be about 25% greater 
than mean network flow for the Poisson line process. 

Of course the second-order behaviours of the flows are rather different: flow at the centre 
of the Poissonian city inherits asymptotically non-degenerate randomness from the random 
configuration of the Poisson line process, while a central limit argument shows that the flows 
at the centres of the two kinds of flow in Manhattan cities are asymptotically deterministic. 



5 Complements and conclusion 

In conclusion we present some notes about complements and issues for further research, 
illustrating the potentially rich theory concerning the Poissonian city. 
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Empirical comparisons 



Clearly the Poissonian city does not accurately represent real cities; there will be variation 
both of geometry and of traffic flow. It would be interesting to make empirical comparisons 
with actual street-map and traffic-flow data, both in terms of network distance statistics 
compared with the results of Section 11.21 and (much more demanding from a data-collection 
point of view) in terms of flow statistics compared with the results of Section [3J One would 
expect qualitative agreement at best, rather than quantitative, in view of the strong stochastic 
assumptions implicit in the Poissonian city. Note however that the results on the asymptotic 
statistics of flow at the centre (Theorem [7]) reflect variation across a sample of different cities, 
rather than within a particular city. 



Lower bounds on path length 



Compare the lower bound of Section 12. 3[ Theorem HI with the more general lower bound of 
Ald ous and Kendalll (120081 . Theorem 2), which holds for all connecting networks using total 
network length proportional to n based on patterns of nodes in a square [0, y/n\ 2 satisfying a 
cert ain quantitative equidistribu tion condition (related to an intuitive coupling construction). 



The lAldous and Kendalll (120081 ) result provides an f2(yTogn) bound on excess, whereas The- 
orem H] uses detailed properties of Poisson line process networks to establish a const, x logn 
lower bound. A natural question is, whether there are any network constructions which 
provide sub-logarithmic mean excess for appropriately equidistributed patterns of nodes, or 
whether on the other hand the general lower bound can be improved. 



Analytical characterization of limit 

The stochastic geometric construction of the limit distribution for flow in the centre of a 
Poissonian city using an improper stationary anisotropic Poisson line process (Section 13. 3j) is 
explicit and lends itself to simulation; however it would be helpful also to have an analytic 
expression or at least characterization of the limiting distribution. This seems difficult. Note 
that we can produce a stochastic representation of the moment generating function M(p) 
of the limit distribution in terms of the following probability: consider a Poisson process of 
intensity a of pairs of points on [—1, 1] x (0, oo): then M(— a) is the probability that no pairs 
produced by this process are separated from the origin by lines of the improper stationary 
anisotropic Poisson line process. Simila r representation s are of use in perfect simulation 



of area-interaction point p rocess models ([Kendalll . 119971 ) and exact simulation of diffusions 



(IBeskos and Roberts! . 120051 ) . However in the current case it is not yet clear whether this offers 
any progress towards simulation methods or delivering a useful analytical representation. 

A slightly easier question is whether the convergence of T n /n 3 to the limit distribution 
holds for all moments. Again, at present no progress on this can be offered beyond the work 
noted in Section 13.41 

Aggregation issues 

What can we say about similar situations where the distribution of nodes generating the flow 
is non-uniform? or even when the nodes generating the flows lie along the Poisson line process 
itself? (Thus precluding the need for the "cross-country" plumbing otherwise requi red to get 



onto the network.) Considerations of this kind are latent in the early work of [Davidson 
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( 119741 ). and it would be interesting to see them applied in the more quantitative set t ing o f 
the present work. It is possible that the coupling construction in lAldous and Kendall! (120081 ) 
would be of use here. 



Three dimensions and higher 

In higher dimensions one needs to consider what kind of network is being deployed. One might 
for example consider the edge process of a Poisson hyperplane tessellation, or alternatively 
one might consider geodesies constrained to lie on the union of all faces of the tessellation. In 
the second case one can derive upper bounds on excess by considering the derivative planar 
problem obtained by taking a 2-plane slice through the source and destination nodes; it is 
then a question how much the excess may be reduced by varying the orientation of the slice, 
and it is a further question whether the excess can be further substantially reduced by using 
paths which do not lie wholly on the slicing plane. It may be possible to make progress in 
the first case by adopting the growth process approach of Section 12.21 



Note that iBoroczky and Schneider! (120081 ) describe higher-dimensional results for similar 
problems; however their results concern standard stereological quantities, while we would 
need results involving infima of lengths of regular curves on the boundaries of Poisson cells. 



Moving beyond line processes 

Certainly one can conceive of results for situations based on processes which approximate 
Poisson line processes; Boolean models based on long line segments, or fibre processes for 
which there is strong control of total fibre curvature. It would be particularly interesting to 
determine the extent to which Poissonian cities and Manhattan cities represent two extremes 
of a suitable class of models. 



User equilibrium 



The notion of UE (User-Equilibrium IWardropl . Il952l . contemporary with the related notion 
of Nash equilibrium) supposes that each user has a utility structure for choosing which 
route they might take based on travel time, which is affected only on available route-lengths 
but also on flow along the routes. Interest is then focused on systems of choices by users 
which result in User Equilibrium; no one user can obtain a shorter route by varying their 
own route. Explorations have already been mad e in the context of queueing theory: see 
for example ICalvert. Solomom and Ziedind (119971). who consider the effect of augmenting a 
simple queueing network and lAfimeimounga. Solomon, and Ziedinsl ( 120051 ) , who consider a 
system of interactions between a -/M/l queue and a -/N^/oo batch queue. There are 
interesting possibilities in the context of the Poissonian city, for example considering that 
traffic from p~ to p + chooses each of the two possible routes prescribed by the semi-perimeter 
algorithm according to considerations both of length and of integrated total flow along the 
routes. 

Such problems are naturally formulated in terms of phase transitions in statistical me- 
chanics, perhaps using the improper stationary anisotropic Poisson line process of Section 
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